


































SUMMARY 



^ The accuracy of sedimentary interval velocities 
derived from a towed hydrophone array can be predicted in the 
ideal case of horizontal interfaces and an idealized source 
pulse. The parameters of the prediction model are the two-way 
normal-incidence travel time and RMS velocity to the horizon of 
interest, array length, and source pulse width. First, the 
range of RMS velocities likely to be encountered in the upper 
500 m of sediments is determined. Second, the error in these 
RMS velocity measurements is calculated as a function of array 
length, height above the sea floor and pulse width. Third, 
the RMS velocity error is converted to interval velocity error 
as a function of layer thickness and velocity, pulse width, 
and array length and height. It should be emphasized that these 
errors do not include any additional errors generated by the 
hyperbolic travel-time assumption and do not include the degrada¬ 
tion encountered for dipping layers. In addition, our discussion 
refers to the analysis of one Common Shot Point or Common 
Depth Point Gather. Considerable improvement in an average 
interval velocity determination can be obtained by the analysis 
of several nearby CSP's or CDP's. This improvement will depend 
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on lateral continuity for the horizons of interest. In the 
discussion which follows, we emphasize that only one CDP or 
CSP is being considered. 


A short pulse, long array, small height above the 
sea floor all will contribute to increased interval velocity 
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resolution. A 1 km long array and a source with a 2 millisecond 
pulse width, both towed 100 meters above the sea floor, can yield 
interval velocities with errors of 2 percent for beds extending 
from 50 to 500 m below the sea floor. For beds thinner than 
50 meters the error increases rapidly, becoming approximately 
10 percent for 20 ra beds. For beds greater than 200 m thick 
whose tops are 400 meters below the sea floor the error is 8 
percent. Again the error increases rapidly for thinner beds. 

For an array altitude of 500 m the minimum error is 7 percent 
and this increases rapidly for beds thinner than 150 m. 

A 1 km long array, 2 millisecond source pulse, and 
a tow altitude of 100 m is recommended to approach the desired 
1 percent error criterion on interval velocity. These require¬ 
ments are within current technological and ships operational 
capabilities. A deep towed array longer than 1 km might prove 
difficult to deploy, tow accurately and safely, and retrieve. 

A pulse width of 2 milliseconds is a compromise between the 
opposing requirements of velocity resolution and depth of 
penetration. A tow altitude of 100 rn might be maintained in 
areas of known topography without the risk of dragging the source 
and array along the bottom. Operational considerations might 
alter the above array and source specification, but the resultant 
change in interval velocity accuracy can be predicted. 

The 1 km array should consist of 25 hydrophones (or 
hydrophone groups) spaced 40 meters apart. The sampling interval 
required is 0.5 millisecond per hydrophone with 12 bit accuracy. 
This yields a data rate of 600 kilobits/second (exclusive of 
engineering data). The relative altitude of the individual 
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hydrophones with respect to the source should be known to 1 m 
and deviations of the array from the horizontal (kiting) should 
be minimized. Because CDP gathering of the data is recommended 
the altitude of the array should not change rapidly over a 
period of 15 minutes. This source and receiver array towed near 
the bottom will yield considerably more accurate interval 
velocities in thinner beds than those obtainable with a surface 
source and array. 
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1. Array Velocities 

Subsurface interval velocities are typically derived 
from array velocity measurements. Thus, the accuracy of 
interval velocity measurements are dependent on the accuracy 
of the determination of array velocities and event arrival 
times. To determine the proper array configuration for a 
deep towed source and receiving array, when the objective is 
to determine accurate interval velocities to a depth of 500 
meters beneath the sea floor, we have investigated the RMS 
velocities expected for various array altitudes above the 
sea floor. (Array velocities are usually equated with RMS 
velocities even though this is only the case when there is no 
dip.) To generalize the analysis we have used an interval 
velocity function which varies linearly with depth. That is 

v » a + bz. 

Since RMS velocities are dependent on the actual interval 
velocities and hence the geology, the use of a linear interval 
velocity depth function where the slope, b, varies from 0 to 
2.5 km/sec/ktn, will include most sedimentary cases of interest. 

In Figure 1, we display the array velocities derived by integrating 
this linear interval velocity depth function. In Figure la, 
for example, the array is at an altitude of 100 m above the sea 
floor, and array velocities were calculated for linear interval 
velocity functions where the slope increased from 0 to 2.5 
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km/sec/km in increments of .1 km/sec/km. The asterisks on 
the array velocity functions indicate depths of 500, 1,000, 

1,500, 2,000 and 2,500 m respectively. At the recommended 
altitude of 100 m the array velocity varies from 1,500 to 3,000 
m/sec, indicating that a broad suite of array velocities are 
available from which interval velocities may be derived. (The 
upper limit is 2,000 m/sec if we consider only the upper 500 m 
of sediments.) 

In Figure lb the array altitude was increased to 
250 n and the suite of array velocities available has decreased. 

In Figure lc, which corresponds to an array altitude of 500 m, 
we see that the range of array velocities available has further 
decreased. In fact, at a tow altitude of 500 m all the array 
velocities for a reflector 500 m beneath the sea floor are 
within the region 1.5 to 1.8 km/sec. To compare the type of 
array velocities that would be available from a conventional 
surface source and receiver array, we have continued calculating 
array velocities for an array altitude of 2,000 and 4,000 m 
above the sea floor (see Figures Id and le). In both cases the 
suite of array velocities available decreases 3till further. 

In particular, when the array altitude (water depth for a surface 
array) is 4,000 ta all of the interval velocities for a subsurface 
reflector 500 m below the sea floor have RMS velocities that 
range between 1.5 and 1.6 km/sec. Thus, an interval velocity 
determination based on array velocity discrimination where the 
array is 4,000 m above the sea floor becomes exceedingly difficult. 

Figure 2 is a detailed plot of array velocities 
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available for an array altitude of 100 m above the sea floor. 

In this plot the array velocities vary from 1.5 to approximately 
2.0 km/sec. This will be the array velocity region of interest. 
Thus the problem of determining the optimum array design for a 
deep towed source and array at an altitude of 100 m above the 
sea floor corresponds to determining the resolution possible 
for array velocities from 1.5 to 2.0 km/sec for reflection 
times up to 1.0 second. 
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2. Array Velocity Resolution 

To develop a measure for defining array resolving 
power, we have considered the case of perfect bandwidth. That 
is, all frequencies for a specified paasbnnd are equally excited. 
Since bandwidth is inversely proportional to the pulse duration 
we consider bandwidth as one requirement for determining the 
resolving power of an array for various array velocities. 

Other variables are array length, reflection depth and the 
actual array velocity. 

Since array velocities are usually derived from an 
approximately hyperbolic travel time—distance relation it is 
common to search through the observed reflection travel time 
data on a trial and error basis for all reasonable hyperbolic 
travel time paths across the array. That is, the two way travel 
time, T, to a reflection event is defined as 

T z = T 0 Z + xVv x 

where T c is the normal two way travel time to the reflection 
event, X is the source-receiver offset and V is the RMS velocity 
as defined by Dix 1955, or in the case of dip as defined by Shah 
1973. By scanning the observed array data for coherent arrivals 
at all the possible two way normal times and RMS velocities it 
is possible to derive ar. interval velocity function using just 
the coherent array arrivals and Dix's 1955 small angle formula¬ 
tion. Thus, before considering the question of interval velocity 
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resolution, we must first consider array velocity resolution. 

(Ue will assume that there is no dip and relatively small source- 
receiver offset so that the RMS and array velocities can be 
equated. Realizing, of course, that if both these conditions 
are not satisfied any interval velocity determination based on 
these assumptions will be in error.) 

In the case of perfect bandwidth where the corresponding 
time resolution is one sample, a direct measure of array resolv¬ 
ing power is the unit sample semblance statistic. Semblance, 
a widely used coherency statistic, is commonly used to determine 
array velocities and arrival times. Usually, semblance is defined 
as the sum of all possible cross correlations between the seismic 
traces for a trial x-t trajectory normalized by the auto-correla¬ 
tions of each seismic trace. Thus, the definition of semblance is 

s = i 

w 4 w t 

where the X- are the data samples across the array for a trial 
trajectory and U is the number of seismic traces. In this 
definition the summation over the correlation window, w, placed 
outside both the upper and lower 3u:na is used to increase the 
energy and therefore the statistical stability. To define array 
resolving power wc will not include these additional summations. 
Rather, we consider that the observed data consists of perfect 
delta functions in time and that array resolving power can be 
based on a unit sample semblance coefficient above a certain 
threshold. This unit sample semblance is a particularly 
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interesting measure of coherency since it is normalized between 
zero and unity (unity being the case of perfect coherence) and 
can be related directly to the variance of the distribution. 
That is, semblance is equal to 


s 



rj-r 

ytTl 


where is the variance and is the second moment of the 

distribution. Therefore, knowing the value of semblance we can 
in fact infer the variance and hence standard deviation of the 
distribution. 

l. T e consider array velocity resolution for the same two 
way travel time, To, as the ability to discriminate between 
array velocities above and below the true velocity. We note 
that for velocities close to the true velocity, many of the small 
offset arrivals will be time coincident for sampled data. Only 
at some offset, x, will the hyperbolic travel time paths diverge 
from the true path. At this point, these time samples will no 
longer contribute to the unit sample semblance statistic and 
we will begin to discriminate between these array velocities. 

For analysis purposes we have used the semblance value of 0.5 
as a measure of discrimination. This implies that we can 
discriminate between coefficients of 1 and 0.5. Using this 
definition the problem of defining array velocity resolution 
is reduced to determining when half of the reflection arrivals 
fall off of a trial hyperbolic trajectory. Clearly this is an 
idealized case and degradation due to additive noise or non- 
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perfect bandwidth will decrease the array resolving power. 

However, this definition will serve as a guide to determining 
the proper array configuration for a given bandwidth. For a 
fixed array length R>!S velocity resolution could of course be 
improved by a large offset between the source and the array. 
However, the accuracy of the vertical incidence arrival time 
would decrease. Since interval velocity is calculated from 
vertical incidence arrival times as well as PUS velocities, it 
is unlikely that interval velocity resolution would improve by 
offsetting the source and the array. 

Our decision to use a semblance value of 0.5 as a 
measure of discrimination is arbitrary, but suitable for the 
following analysis. The actual array length required will then 
be equal to twice the offset at which the arrivals first begin 
to fall off the trial trajectory. (Any other semblance value 
would imply an array length of 1/S of this value.) Semblance 
as a measure of co crency is also sensitive to the number of 
estimates available. For example, in the noise free case if only 
two estimates are available (two channels), the semblance statistic 
is a poor discriminator since the only possible values are 0, 

.5, and 1. As the number of channels increases the discrimination 
capability of semblance improves. (Where noise alone is present, 
no coherent arrivals, we expect the summations in the numerator 
and denominator to be equal and the limiting semblance value 
would be 1/H.) Thus, if a higher semblance value is used for 
a basis of discrimination, tiie number of channels that would 
be required should also increase to maintain the comparable 
level of discrimination. 
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To determine the offset v?hcre the arrivals begin to 
diverge in the case of perfect bandwidth, we take the partial 
derivative of time with respect to array velocity for the 
standard hyperbolic travel time equation and solve this equation 
for the array size necessary for a specified percent error, £, 
in array velocity. That is, 

x ~ 0 e j ( 

v:here i’ is the two way normal time, B is the bandwidth and 
V is the array velocity. In Figures 3, 4, and 5 we have displayed 
the array size necessary for a specified hi! 3 velocity error at 
vertical two way travel times of .133, .333, and 1.333 seconds 
be lew the array for bandwidths of 100, 200, 500 and 1,000 Hz 
based on a semblance coefficient of 0.5. (In these plots the 
array length calculated from the above formula was doubled to 
determine the actual array length required.) In Figure 3a wc 
display the array size necessary for u vertical two way travel 
titie of .133 second and a bandvv’idth of 100. To achieve a one 
percent error in array velocity determination, an array size of 
3.9 Un is required for an array velocity of 2.0 km/sec. In 
Figure 3b the bandwidth increases and an array size of about 2.1 
km is required to achieve the same percent array velocity resolu¬ 
tion. In Figure 3c, where the bandwidth has now become 500 Hz 
it is possible to achieve the same array resolution with an array 
size of 1.0 km. In Figure 3d where the bandwidth is 1,000 Hz 
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we see that it is possible to achieve the same percent error 
with an array sice of approximately 0.6 km. Tlius, as the band¬ 
width increases the required spatial aperture decreases. Clearly 
it is desirable to tow the minimum array length and to have as 
broad a bandwidth as possible. However, achievement of the 
necessary bandwidth will be restricted by the spectral 
characteristics of available sources. 

In Figures 4 and 5 we perform the same analysis but 

* 

the time to the reflection horizon is increased. As expected, 
when the time to the reflecting horizon beneath the array 
increases the required bandwidth and/or the array size must 
increase to achieve the same accuracy in velocity resolution. 

Based on Figure 2, we expect the array velocities to 
be between 1.5 and 2.0 km/sec for titles of approximately 1 to 
1.5 seconds beneath the array corresponding to depths of 500 m 
below the sea floor for reasonable interval velocity functions. 

In addition, we expect that it is possible to obtain source 
bandwidths on the order of 200 or possibly even 500 Hz. This 
would indicate, based on Figures 3b and c, that an array size 
of 1 km will be adequate for resolving array velocities to within 
one percent for the expected arrival times and for velocities 
associated with sedimentary horizons up to 500 m below the sea 
floor. Clearly, if more bandwidth can be obtained and the array 
size increased, the resolution will also increase. The bandwidths 
we have considered are idealized in the sense that we assume that 
all frequencies are equally excited within this band. A bandwidth 
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3. Spatial Aliasing 

Spatial aliasing is encountered when the sampling 
interval in space is too coarse for a given phase velocity 
across the array. For example, in Figure 6 we plot frequency 
versus wave number and have indicated a phase velocity of 
2 km/sec. We have also indicated the Hyquist wave numbers 
corresponding to receiver separations of 100, 50, 33 and 25 m. 
Whenever energy is traversing the array at a phase velocity of 
2 km/sec it will be aliased and appear to lie along the lines 
slanting upward to the left. For instance, at 25 m spacing 
a 2 km/sec phase velocity is aliased above 40 Hz. This would be 
a problem if frequency-wave number processing was to be performed 
on the original array data. Usually this is not the case. It 
is common practice to scan the array arrivals only in a specified 
array velocity band since one knows approximately the array 
velocities to be determined. Thus, even though the data may be 
aliased it is often unimportant in practice. If one were to 
process the original array data in the frequency-wave number 
domain (and were to discriminate based on phase velocity) the 
alias lines would have to be followed in the manner indicated 
in Figure 6. 

Dense spatial sampling is required to remove aliasing. 
In Figure 6 a spatial sampling interval of 10 m is required to 
have unaliased data above 100 Hz. Since this necessitates 
increasing the number of channels, it becomes difficult to 
achieve this type of spatial resolution. The exploration 
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industry is rapidly moving towards longer arrays with denser 
spatial sampling and this has necessitated recording on such media 
as video tape because of the correspondingly high data rates. 

Day to day exploration, however, is commonly carried out with 
100 or 50 m array element spacing intervals. While higher 
spatial resolution is clearly desirable, this v.’ould increase 
the number of channels and require significantly greater digiti¬ 
zation bandwidth than is currently available. 






"1 

i 


-16- 


] 

"! 

J 

j 

] 


] 


1 

J 

—» 

J 


4. CDP Versus Common Shot Data 


Array data is acquired in a Common Shot Point (CSP) 
mode, that is, all channels are recorded for a given shot. In 
exploration, this data is reordered into the Common Reflection 
Point or Common Depth Point (CDP) mode. CGP and CDP are 
equivalent geometries in the case of no dip. In the presence 
of dip, however, the CDP geometry has significant advantages. 
Basically this geometry averages the ray paths such that the 
PC'S velocity determined from the array velocity can be used 
to give a better indication of the subsurface interval velocities. 
Several papers have been written, for example, Shah 1973, which 
indicate how array velocity determinations can be turned into 
true RMS velocity determinations even in the presence of dip 
for the CDP geometry. Thus, by measuring the time dip on a 
normal incidence record section and using array velocities one 
can correct for the presence of dip to improve the interval 
velocity determination. In addition, in most cases the effect 
of modest dips on the CDP geometry are quite minor. CDP 
geometry has a significant advantage over CSP geometry for 
horizontally discontinuous or rough reflectors. In the CSP 
mode an interface would have to persist laterally for at least 
one half the array length in order that true reflections from 
it be recorded on all hydrophone groups. In the CD? mode this 
persistance requirement is reduced to the smear in the CDP 
caused by errors in the Shot Point placement (navigation) and 
reflection dip. This advantage becomes important when attempting 
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to measure RMS velocities from rough interfaces such as oceanic 
basement, since one need only find a short segment where basement 
is smooth and horizontal. Additionally, after applying moveout 
corrections and stacking the section will be a better representa¬ 
tion of the sub-sea-floor geology in the CDP case. 

To acquire CDP data using a 1 km deep-towed array 
with 25 channels and 40 m spacing would require a 20 m shot 
spacing interval for full 25 fold CDP coverage. If the array 
were towed at 1 kt (.5 tn/sec) then a 40 sec repetition rate 
would be required. It would be preferable to tow the array at 
2 kt or 1 m/sec 3nd in this case a 20 sec repetition rate would 
be i-equired. The CDP mode requires accurate positioning of the 
shot points to avoid smearing the Common Reflection Points 
(CRP). For 25 fold data, 25 consecutive shots contribute to 
each CDP. It is the average spacing of the shot3 that must be 
controlled accurately, rather than the spacing between individual 
shots. For a 1 km array the 25 shots are spaced over a distance 
of 430 m. If one specifies a maximum smear in the CRP's, 
caused by misplacement of the shot points, of 5 m,- then the 
average shot spacing must be accurate to 1 percent. It is 
assumed that longitudinal deformation or stretching of the array 
will be small (i.c, 'v2 m). If the array is decoupled from the 
source, variations in the distance between the source and the 
array must be minimized and known for each shot point. For data 
acquired at 1 kt it would take 1,000 seconds or 16 minutes to 
acquire one CDP. Therefore, only slow deformations of the array 
would be tolerable. At 2 kts, however, the data would be 
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acquired in 8 minutes and more rapid deformations could be 
tolerated. 

To acquire CDP data with a deep towed array will 
require that kiting be kept to a minimum. The minimum 
requirements of 500 Hz bandwidth necessitates that time be 
known to less than a millisecond, that is, we must know the 
array height to within a meter. In one meter a pressure 
change of 1.5 psi occurs and although an absolute measurement of 
pressure is not necessary, relative pressure and depth should 
be recorded to the required accuracy. (An accurate calibration 
of the pressure sensing units would be necessary.) 

Two corrections are required because of kiting. Hie 
first is a timing correction and the second is a spatial 
correction. The array appears smaller as the kiting angle 
increases, thus it is absolutely necessary that the array 
deformation be known and that this deformation be removed on 
a shot basis prior to the CDP gather. Instrumenting the array 
with pressure sensors that have been calibrated initially 
and have an accuracy to better than 1.5 psi should provide 
sufficient information to remove the effect of kiting. In 
addition, if the array is towed with a drogue and the array 
is neutrally bouyant and the source and array are mechanically 
isolated from the towing cable, the corrections after the initial 
settling should be a minimum. 
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5. Interval Velocity Error Estimation 


In a previous section we discussed array velocity 
errors. For relatively snail offsets and no dip they can 
be considered RMS velocity errors. Here we relate RMS velocity 
errors to interval velocity errors. 

Interval velocities are calculated from pairs of RMS 
velocities and their associated vertical incidence two-way 
travel times via Dix's 1955 formula: « 


v: 
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v* r v - v, l r> 
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where and are respectively the RMS velocity and travel 
time to the bottom of the layer, and and refer to the 
top of the layer. Each of these terms has an error associated 
with it. We know how to calculate RMS velocity errors and we 
assume that travel-time to a reflector can be determined to 
an accuracy of one source pulse width (i.e. the reciprocal of 
the bandwidth). From known errors in V^, T^, and there 
are several ways to estimate the error in'V',... 

IN 1 

One could assume that the errors in each of the terms 
add algebraically. In this case we would add percent errors 
during a multiplication or division and add actual error for 
additions or subtractions. The square root requires halving 
the percent error. This i3 a form of worst-case error 
estimation, since we assume that each term is in error by the 
maximum amount and in the most harmful direction. 
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A second method assumes that the errors add vectorially. 
Instead of adding errors (either percent or actual) we assume 
that they are at right angles to each other. Thus an error of 
3 percent and an error of 4 percent would yield a combined error 
of 5 percent (3^ + 4^ = 5^). This is a form of most probable 
error, since we assume that each term will not always be in error 
by its maximum amount, that is the errors are independent of each 
other. 

We have chosen a third method 'which is also a form of 

worst-case error estimation. First we calculate by assuming 

that is larger than and is smaller than their true values 

by the maximum allowable error. Thus if V is 2,000 m/sec and 

the error estimate is 1 percent we set equal to 2020 m/sec. 

Similarly if T is 1.0 second and the error is + 2 milliseconds 
2 

(1/500 liz) then T^ is set to 0.998 seconds. In a similar manner 

we underestimate and overestimate T^. This will yield a large 

value of V w ..,, see Figure 7. We then calculate low values for 

(1980 m/sec) and T^ and high values for and (0.993 sec). 

This gives us a low value of V TKI . Then the percent error in V 

MIN 

is estimated as: 


r - 'V -V~- 

to*. ^ ,v 


5-^-r 


With this type of worst-case error estimation, the actual errors 
will be less than the estimated errors. However, wc have pre¬ 
scribed error limits and the true interval velocity will be 
with.in these limits. 
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6. Interval Velocity Resolution 
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The above discussion gives a proper background to 
interpreting Figure 8. In Figure 8 we plot interval velocity 
percent error as a function of layer thickness (Delta Z). (It 
must be emphasized at this point that these errors do not include 
any error associated with assuming hyperbolic trajectories.) 

Figures Sa to e show the interval velocity error as a 
function of layer thickness for a 1 km array towed 100, 250 and 
500 m above the layer for bandwidths of 500 and 200 Hz. The 
six curves on each plot are for layer velocities from 1,500 
m/sec to 2,000 m/sec in 100 m/sec steps. The material above 
the layer is assumed to have a velocity of 1,500 m/sec. The most 
important feature of these plots is the strong dependance of the 
error on array altitude and bandwidth. These emphasize the need 
for a large bandwidth source and a low towing altitude. Low 
altitude not only achieves higher velocity resolution for thick 
beds, but extends this high resolution to thin beds (down to 100 m 
in Figure 8a). It appears that 2 percent accuracy can be 
achieved with a bandwidth of 500 Hz and an altitude of 100 m. 

Two other features are evident; the relatively small effect 
interval velocity has on resoition, and in Figures 8a to 8d 
a slight decrease in resolution for thicker beds. This decrease 
in resolution results from flatter hyperbolic trajectories at 
the base of the bed as the time to the base increases. This 
flattening increases the error in the RMS velocity. 
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By way of comparison with a surface source and array 
we have included similar plots for arrays of 4.8 Ion and 10 km 
length in a water depth of 4 km. (Figures 8f and g) The 10 km 
array is hypothetical as we do not know of any in existence. 

Note that these two plots have a different thickness axis. For 
a 1 km thick bed the 4.8 km array has a minimum error of 3 percent 
and a 10 km array improves the accuracy to 1 percent. However, 
the velocity error increases rapidly as the bed thickness decreases 
below 500 m. ' 

In conclusion we see that interval velocity accuracies 
near 2 percent can be obtained with a 1 km array towed 100 m 
above the sea floor with a source of at least 500 Hz bandwidth. 
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7. Refractions 


Critical angle refractions or head waves might be 
observed with a 1 km array if it is towed close to the sea floor 
and there are shallow high velocity beds. These refractions 
would px*ovide a useful measure of velocity at the top of the 
refracting horizon. Refractions normally observed with surface 
sonobuoys are generally low frequency which makes it difficult 
to accurately measure their arrival time across the array. 

The 1 km array coupled with accurate timing could yield good 
phase velocity estimates. However, in doing RMS velocity 
scans it may be necessary to mute or zero the large offset 
traces since these arrivals will lead to a degradation in the 
RMS velocities and thus the interval velocity determination. 
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Figure C ap tio ns 

Figure 1. Array velocities derived fro.n linear interval velocity 
vs. depth functions for array altitudes of 100, 250, 500, 2,000, 
and A,000 m. On each plot the normal two way travel times 
corresponding to depths of 500, 1,000, 1,500, 2,000, and 2,500 m 
are indicated. For an array altitude of 100 m the expected array 
velocities range from 1.5 to 2.0 km/sec and the expected two way 
normal times range from 0.5 to 1.0 sea;. Increasing the array 
altitude decreases the suite of array velocities and increases 
the two way normal time. For an array altitude of 4,000 m all 
array velocities for a reflector 500 m below the sea floor are 
between 1.5 and 1.6 km/sec. 

Figure 2. A detailed display of the expected array velocities 
for an array altitude of 100 rr.. 

Figure 3. Array size (km) vs. percent array velocity error for 
a reflector at .133 sec of two way normal time (depth of 100 m 
for a velocity of 1.5 km/sec) for bandwidths of 100, 200, 500 
and 1,000 Jlz and for array velocities 1.5, 2.0, 2.5 and 3 km/sec. 
For a fixed two way normal time, as the bandwidth increases the 
percent array velocity error decreases. 


Figure 4. Array size (km) vs. percent array velocity error for 
a reflector at .333 sec of two way normal time (a depth of 250 in 
for a velocity of 1.5 km/sec) for bandwidths of 100, 200, 500 and 
1,000 Hz and for array velocities 1.5, 2.0, 2.5 and 3 km/sec. 






Figure 5. Array size (krn) vs. percent array velocity error for 
a reflector at 1.333 sec of two way normal time (a depth of 
1000 m for a velocity of 1.5 km/sec) for bandwidths of 100, 

200, 500 and 1,000 Hz and for array velocities 1.5, 2.0, 2.5 
and 3 km/sec. 
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Figure 6. Frequency vs. wave number plot with a phase velocity 
of 2 km/sec indicated. The Nyquist wave numbers corresponding 
to spatial sampling intervals of 25, 33, 50 and 100 m are indicated. 
Ifnenever energy is traversing the array at a phase velocity of 
2 km/sec it will be aliased and appear to lie along the lines 
slanting upward to the left. For the 2 km/scc phase velocity 
and a 25 m spatial sampling interval all frequencies above 40 Hz 
will be aliased. 

Figure 7. Diagram indicating the method used to compute the percent 
interval velocity error from RMb’ velocity and normal two way travel 
time measurements. To determine the minimum interval velocity we 
overestimate the RMS velocity and underestimate the time at the 
top of the layer and underes timate the RMS velocity and overestimate 
the time at the bottom of the layer. To estimate the maximum 
interval velocity we underestimate the RMS velocity and over¬ 
estimate the time at the top of the layer and overestimate the 
RKF velocity and underestimate the time at the bottom of the layer. 
Half the difference between these two interval velocities divided 
by the true interval velocity is the percent error used in Figure 8. 



Fig,uras 8j to c. Percent interval velocity error vs. thickness 
(Delta 2), for an array of 1 km towed 100, 250 and 500 m above 
the sea floor and a source bandwidth of 500 Hz. Each curve 

corresponds to an interval velocity of 1.5, 1.6, 1.7, 1.8, 1.9 and 

2.0 k:n/sec respectively. As the array altitude increases the 
percent interval velocity error also increases. As the thick¬ 
ness of the beds decreases the percent interval velocity error 
increases significantly in all cases. As the thickness of the 

beds increases the percent interval velocity error eventually 

% 

will increase because of the increased array velocity error. 

Figures 8d and e. Percent interval velocity error vs. thickness 
(Delta Z) , for an array length of 1 km towed 100 and 250 rn above 
the sea floor and a source bandwidth of 200 Hz. Each curve 

corresponds to an interval velocity of 1.5, 1.6, 1.7, 1.8, 1.9 and 

2.0 km/sec respectively. In comparison with Figures 8a to c the 
decreased bandwidths result in significantly increased percent 
interval velocity error. 

Figure 8f and g. Percent interval velocity error vs. thickness 
(Delta Z), for an array altitude of 4000 m, source bandwidth 
of 500 Hz and array lengths of 4.8 and 10 km. Each curve corresponds 
to an interval velocity of 1.5, 1.6, 1.7, 1.8, 1.9 and 2.0 km/sec 
respectively. For this array altitude (or for a surface array 
where the water depth is 4000 m) a 3 to 4 percent interval velocity 
error can be obtained for layers with a thickness of about 1 
km for an array length of about 4.8 km. For an array length of 
10 km the percent error decreases significantly. For both array 
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DIMENSION XC 2000 ),ERC 2003) 

C****™****™********'************-******#************************** 

C DETERMINE REQUIRED ARRAY LENGTH FOR GIUEN TARGET DEPTH, 

C BANDWIDTH & RMS UELOCITr 
C ENTER PARAMETERS 
CALL ERASE 

ACC£PT"T0 (MSECS) ? ",T0 
T0-T0X1000 

ACCEPT "URI1S0, URMS MAX, DORMS (fVSEC) URMS0, URMSX, BURNS 
NU-CURMSX-URMS0)/DURMS+i 
EURMS-DURMSX1000 
URMS0-URMS0X1000 

ACCEPT"U X ERROR MIN, MAX ? ", CHIN, EMAX 

ACCEPT"SEM THRESHOLD (E G. .5) ? ", SMIN 

XMPL-4.; MAX ARRAY FOR PLOT 

DELER-.05 

EMIN-EMIN'100 

D-1AX-EMAX/100 

DELER-DELER/100 

ACCEPT"X ERROR ( IN XX E.G. I ), X-AXIS ( INxRM EG. 5)7 ".DX.XSCAL 

NE“(EMAX-EMIN)/DELER+i; # OF ERROR STEPS 

DO 1000 jPLOT-1, X, LOOP FOR 4 BANDWIDTH5 

CALL PLOTS(IDUM, “TAXUF") * 

CALL PLOT( 3. , 3 ., -3) 

IFCJPLOT.EG.1) FCUT“100 
IFCJPLOT.EQ.2) FCUT-200 
IF( JPL0T .E0.3) FCUT-500 
IFCJPL0T.EQ.4) FCUT-1000 

DT-l./FCUT, TIME DURATION OF PULSE 

£+*+***,#****#*##******##**■**********'*****+*»#****#+##*’** ******** 

C DRAW A>£S 

CALL AXISC 0 ,0 , 7HX ERROR, -7, DX*1M?*NE*DELER,0. , 0., 1 /DX) 

XL*=XSCAL*XMPL 

IFCXL.GT.10) XL-XSCAL 

CALL AXISC 0.,0 , SHARPAY, S, XL,90.,0.,1.xXSCAL) 

CALL PLOTC0. ,0. ,3) 

c*+******.*****#****w*****r*#**#**r***#*r*w#**w##$*###*M***#4***i«*# 

C LOOP FOR ALL RMS UELOCITIES 
DO 200 Ji-l.NU 

URMS-URMS3+C JJ-1 )*EVRMS, NEW URMS 

JER-1 

DO 50 J-l,NE; LOOP FOR ALL ERRORS 

EL'-EMIN+C J-l )*DCLER, GET ERROR THIS CASE 

TT-T0*T0xC DT*DT) 

TT-EU*EU*TT 
RADi-TT+0 2S 
RAD1-S0RTCRAD1) 

RAD2-RAD1+0 5 
RALOSOR r C RAD2) 

XCJER)-URMS*DT*SAD2xEU, X FROM EO, NO SEM CONDITION YET 

XCJER)-XCJERVSMIH, USE SEM MIN TO GET REQUIRED ARRAY 3IZ 

ERfJER)-EU*100; X ERROR 

IFCXCJER) GT.XMPL) GO TO 50 

IFCJER EQ 1) EC-CU, SPUE 1ST ERROR ON PLOT 

JEND-JER 
JER-JEP+1 
50 CONT WLE 

C+*<r ******* *«****»» *************** ******************a************* 

C PLOT RESULTS THIS CASE 
IPEN-3 

DO 100 K-l.JEND 





XX-ia0*E0*DX+100*(K-l)*DX*DELER 

YV-XCK)*XSCAL 

CALL PLOT(XX, YY, IFEN) 

100 IFEN-2 

CALL NUMBER(XX+ 2S, YY, . 1, URMS, 0. , 3) 

CALL PLOT(0. ,0. ,3) 

200 CONTINUE 

C CLOSE PLOT FILE, PLOT 

CALL PLOTC0. ,-2. ,3) 

CALL SYMBOL (939.,999. , .2, 4H TO ,0,4) 

CALL NUMBERC999. ,999. , .2, T0, 0 ,2) 

CALL SYMBOLC 993.,999. , .2, 14H BANDWIDTH ,0 ,14) 
CALL NUMBER(999. ,999. , . 2, FCUT, 0. , 1) 

CALL PLOT(0. ,0. ,999) 

CALL FStf'PC 'YX^LT SO" ) 

DO 500 JJJ-1, 1 

CALL FSWAPC "FPBITS. SU") 

CALL FPUOK 
503 CONTINUE 
1003 CONTINUE 

CALL RESET 
CALL EXIT 
END 




Appendix II: Interval Velocity Error Computation 


C DETERMINE INTERVAL VELOCITY ERROR (X) 

C FROM ARRAY VELOCITY & TWO WAY TRAVEL TIME MEASUREMENTS 
C ENTER PARAMETERS' 

CALL ERASE 

ACCEPT"BANDWIDTH (HZ) ? ", BAND 
ACCEPT"APRAY LENGTH ? ", XARAY 
ACCEPT"ARRAY ALTITUDE (KM) ? ", AH 

ACCEPT "INTERVAL: UMIN. UMAX, DELU(E.G. I. 5,2., .1) >? '*, VMN, U1AX D* 
ACCEPT"0-TIME, 1-DEPTH ? ",IDTP 

IFdDTP.EQ.0) ACCEPT "MAXIMUM LAYER THICKNESS (SEC) ? ", TH 
IF'IDTP. EQ. 1) ACCEPT "MAXIMUM LAYER THICKNESS (KM) ? ", TH 
IF(IDTP.EO.a) ACCEPT"INKSEC (EG. S) ? ", DX 
IF( IDTP.EO. 1) ACCEPT"IN/KM (E G. S) ? ", DX 
ACCEPT”MAX X ERROR, INK TERROR ? ", EMAX, DY 
N-TH*EAND, i OF STEPS 

DELT-1.xBAND, DELTA T 

NL'-CVMAX-VMNVDV+l; ♦ OF VELOCITIES 

C^'*#***».*****#*#*#^**W#****#**#*it-i(«-’i(»**Y*)t(#<»-*****)t:**>'Ji!t(><» t )|(*** : i((» 

C OPEN PLOT FILES. DRAW AXES 

CALL PLOTS(IDUM,“TAXVF” ) 

CALL PLOTO. , 10. ,-3) 

IF( IDTP. EQ. 0 )CALL AXISC0. , 0. , 14HDELTA TO (SEC),-14, TH*DX, 0. , 0. , I -'D> 
IF(IDTP.EQ.l) CALL AXIS(0., 0.,12HDELTA 2 (KM),-12, THKDX,0. , 0. , I.kbx 
CALL AXISC0 , 0. , 7HX ERROR, 7, EMAX*DY, 90. , 0. , 1. kDY) 
C*M****#****M**m*#**#m*m'#M*W*m******************* 

C 1ST LAYER IS FIXED, GET VRMS (VI) & T0 (Tl), URMS ERROR (VER1) 

C & VI LOW (V1L), & VI HIGH (V1H1 
C DEFINE AFPECTIVE ARRAY LENGTH**2 (XAR2) 

1)1*15 

T1-2*AHKV1 

XAR2=( SrXARAY)**!, (ASS'JMSE SEM -.5) 

V12-U1**2 

VERl-DELTW12*SORT( T1**2+XAR2/V12 )KXAR2 
V1L-V12*(1.-UER1 )**2*( Tl+DELT) 

U1H-V12*<1.+VER1)#*2*( Ti-DELT) 

C# ***** ***> **v:*****#***#*************:********-***X******** ******* 

C LOOP FOR ALL VELOCITIES 
DO 100 JV-l.NV 

UINT-VMN+<JV-1)*DV, NEW INTERVAL VELOCITY 

VINT2-VINT**2 

C*#»»»**** i l(>M)*^>M(*»t'**'****il’»G|!'*-**»*.**%*iH»f***«* ****** ****** 

IPEN-3 

IF(IDTP.EQ.I) N“PAND*2*ThUVINT+l 


C****#*^**************^-************************** 


NEXT LAYER THICKNESS 


LOOP OVER ALL THICKNESSES 
S0 DO 90 J-4,N 
DT«J*DELT, 

TZ-T1+DT 

URMS2-(VINT2*DT+V12*T1)KT2; 2ND VRMS (TRUE) 

VER2-DELT*VRMS2*S0RT( T2**2+XAR2xVRMS2)KXAR2; V. ERROR THIS VRMS 


V2H=VRMS2*<1.+VER2>**2*( T2-DELT), 
VMAX"(U2H-V1L)K(T2-Ti-2*DELT), 
V2L-VRM52*<1.-VER2)**2*'(T2+DELT), 


V2 RMS HIGH*#2 V TO - DELTA 
VINT MAX **Z 

V2 RMS LON**2 # 70 + DELTA 


VMIN-( U2L-V1H )K( T2-T1+2*DELT), 
IF(UMAX.LT 0) GO TO 93 
IIXUMIN LT 0) GO TO S3 
L'MAX-SQRT(VMAX) 
VMIN-SGRT(VMIN) 

DJTT-VMAX-VMIN, 

DV1NT-.5TDUTTKUINT, 

1F(DUINT LT.0 ) GO TO 90 


VIITT MIN **2 


DELTA VINT 
X VINT ERROP 


If(DUINT. GT. .0i*EMAX) GO TO 90 

C***V*«***fc*V4ntt*)|&)M¥************#*******»**¥i***********>IOIit**« 

C PLOT * LINT ERROR 

xx«dt*dx 

1F(IDTP EQ. 1) XX-D/*DT*OlNT/2 

yv-i00*duint*dy 

CALL PLOT(XX, YY, IPEN) 

IPEN-2 

90 CONTINUE 

IF(JU.EO.1.OR.JU. CO . NU> CALL NUMBERCXX+. 25, YY-. 05, . 1, UINT, 0 
100 CONTINUE 

C4**:Y44*Y4**4**'***»*****^**^*4**4*'4444*'44» , *iW*4**444.*4*''*4w , it!4*44*'!M( 

C. ANNOTATION, CLOSE PLOTFILE, PLOT 

CALL PLOTC TH4DX/2,EMAXWDY,-3) 

CALL SYM3OL(0. ,0. , .1,15HBANDWIDTH (HZ) ,0. ,15) 

CALL NUMBER(999., 999., .1, BPjffi, 0.,1) 

CALL PLOTC0.,-.2,-3) 

CALL SYM3OLC0. , 0. , .1, 19HARRAY ALTITUDE (M) ,0. ,19) 

CALL NUM3ER(S39.,999.,.1,AH*1000,0 ,1) 

CALL FLOT(0. ,-.2,-3) 

CALL SYMBOL(0. ,0. , .1,18HARRAY LENGTH (KM) ,0. , 18) 

CALL NUMEER(999.,999. , .1, XARAY, 0 ,1) 

CALL FLOTC0. ,0. ,999) 

CALL FSi>JAP( "AGPLT. SV” ) 

CALL FSWAP( "FPBITS. SO'* ) 

CALL EXIT 
END 





